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ABSTRACT 



The picture of the early evolution of globular clusters has been significantly revised in recent years. Current scenarios require at least 
two generations of stars of which the first generation (1G), and therefore also the protocluster cloud, has been much more massive 
than the currently predominating second generation (2G). Fast gas expulsion is thought to unbind the majority of the 1G stars. Gas 
expulsion is also mandatory to remove metal-enriched supernova ejecta, which are not found in the 2G stars. It has long been thought 
that the supernovae themselves are the agent of the gas expulsion, based on crude energetics arguments. Here, we assume that gas 
expulsion happens via the formation of a superbubble, and describe the kinematics by a thin-shell model. We find that supernova- 
driven shells are destroyed by the Rayleigh-Taylor instability before they reach escape speed for all but perhaps the least massive and 
most extended clusters. More power is required to expel the gas, which might plausibly be provided by a coherent onset of accretion 
onto the stellar remnants. The resulting kpc-sized bubbles might be observable in Faraday rotation maps with the planned Square 
Kilometre Array radio telescope against polarised background radio lobes if a globular cluster would happen to form in front of such 
a radio lobe. 
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1. Introduction 

Galactic globular clusters (GCs) today typically consist of old 
low-mass stars and little or no gas. However, there must have 
been a time when they formed as gas-rich objects with numer- 
ous formation of also massive young stars. Many details of this 
early epoch have only recently been discovered dGratton et al.l 
l2012tlCharbonnell 20101 for recent reviews). Progress has in par- 
ticular been made via spectroscopy and chemical and dynamical 
evolution modelling. This has led to a picture of star formation 
in mu ltiple episodes: In summary (e.g. iPrantzos & Charbonnell 
120061 and references therein), the stars in individual GCs are 
mono-metallic regarding the iron group elements (Fe, Ni, Cu), 
and have little scatter and similar trends as field stars for the 
neutron capture (Ba, La, Eu) and the alpha-elements (Si, Ca). 
However, light elements present strong variations from star to 
star with anti-correlations, between O and Na, and Mg and Al, 
respectively. The interpretation is that GCs form from uniformly 
pre-enriched gas, which explains the similarities for the iron 
group, neutron capture, and alpha-elements. To explain the anti- 
correlations, one requires processed mate rial that has been sub - 
ject to hydrogen burning at about 75 MK (IPrantzos et al.ll2007h . 
These conditions are found in the most massive fast rotating 
massive stars (FRMS) and in massive asymptotic giant branch 
(AGB) stars. Thus, one requires a first generation (1G) of stars 
including massive stars, the ejecta of which form a second gen- 
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eration (2G) containing low-mass stars (typically the majority 
of the stars we observe today). The stellar ejecta have to be 
mixed to a varying degree of about 30-50 % with pristine gas 
to produce the abundance patterns (anti-correlation) of the now 
observed 2G stars, but the inclusion of processed gas ejected 
in supernovae (SNe) has to be avoided (|Prantzos et all l2007t 
iDecressin et al 1 l2007bl: ID 7 Ercole et al.ll2.01 ll) . Gas expulsion by 
these SNe seemed to be an obvious way to remove their ejecta 
from the GC together with the bulk of the gas. Both the AGB 
and the FRMS scenario agree that the first generation of stars, 
and thus the initial total stellar population, was much more mas- 
sive than the second generation. If the initial mass function 
(IMF) would have been normal, many of the 1G stars would 
then have had to be lost (becressin et al.ll2007aL IVesperini et aD 
120101: ISchaerer & Charbonnell 1201 lb . Assuming mass segrega- 
tion and formation of the 2G stars in the vicinity of the more 
tightly bound massive 1G stars, a quick change of the gravita- 
tional potential may unbind the major part of the 1G low mass 
stars in the outskirts of a GC. 

Winds and SNe produce interstellar bubb les, as commonly 
obser ved in the interstellar medium (e.g. IChurchwell et al.l 
l2006h . Given the small separations in GCs, they should soon 
unite and thus for m a superbubble (e.g. iBagetakos et ail 1201 It 
iJaskot et"aTl 1201 1|) . GCs are extremely tightly bound systems 
with half-ma ss radii of typically a few, sometimes only one pc 
dHarrislll996[). They are unlikely to have been less concentrated 
in the past ( Wilkinson et alj2003l) . Gravity poses a profound ob- 
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stacle to escaping superbubbles. Superbubbles first need to build 
up a high pressure to lift up the gas. Once the half-mass radius is 
reached, gravity declines quickly, and the pressure force strongly 
dominates, which leads to acceleration of the shell, and thus trig- 
gers the Rayleigh-Taylor (RT) instability. When RT modes of 
about the bubble size are able to grow, the shell fragments and 
releases its internal pressure. This will favour its fall back into 
the central part of the cluster. Here we show that this process 
prevents gas expulsion by SN feedback in all but the least mas- 
sive GCs. The power released by accretion onto dark remnants 
could be sufficient to expel the gas. 

2. Superbubble formation and gas expulsion 

Once an amount of energy comparable to the binding energy is 
liberated, we expect a superbubble to form. The evolution of GC 
superbubbles has been modelled by iBrown. Burkert. & Truranl 
dl99ll 11995 -TT. They showed that the thin shell approximation 
(compare below) models the superbubble expansion faithfully. 

2. 1 . The thin-shell model 

We model the superbubble with the spherically symmetric thin- 
shell approximation, where the change of the shell's momentum 
is simply given by the applied forces, 



d_ 

dt 



(M\>) =pA- Mg . 



(1) 



Here, Ai = 4n PgfV'V 2 d/ is the mass in the shell, with the 
gas density p 2 and the shell radius r; v is the shell velocity, p the 
bubble pressure, assumed to dominate over the ambient pressure, 
A = Anr 2 the surface area of the shell and g the gravitational ac- 
celeration. The bubble pressure is p-(y— l)(rjE(t)- Mv 2 /2)/V, 
with the bubble volume V = 4nP/3, the energy injection law 
E{t), an efficiency parameter 77, and the ratio of specific heats, 
7=5/3. 

As input to the model, we need to specify the mass pro- 
file and the energy input. Following other recent work (e.g. 
iBaumgardt et al.l2008HDecressin et al we use a Plummer 

model for the spatial distribution of gas and stars. The gas mass 
inside a radius r is then given by 



M(r) = (l-e sf ) 



M tot r 3 



(r 2 + r c 2 ) 3 / 2 



where r c = (2 2 ^ 3 - I) 1 ^ 2 ri/2 ~ n/2/1-3 is the core radius and e s f 
the star formation efficiency. For the gravitational acceleration, 
we take into account the stars and half of the gas mass in the 
shell. This results in 



(2) 



1 + 6 sf GM 
l-e sf 2r 2 



(3) 



Our standard case is a protocluster of M tot = 9 x 10 6 M Q , a 
half-mass radius of n/2 = 3 pc, a star formation efficiency of 
1/3, and a Salpeter IMF, which s hould be applicable f or more 
massive GCs such as NGC 6752 dDecressin et al.ll2010l) . 

We have tested three scenarios for the energy injection 
law E(t): In our standard scenario, we assume that all stars 



1 While the 2G formation scenario in the supershell proposed in 
these papers did not stand up to observational scrutiny (e.g. because 
SN ejecta are now thought not to be mixed with the gas that forms the 
2G stars), the hydrodynamics is still valid. 



with initia l masses betwe e n 9 a nd 120 M Q explode as SNe. 
Following iDecressin et al.l (120 1 Ot) . we assume the SNe to con- 
tribute 10 51 ergs, each, with an efficiency of 77 = 0.2. We take into 
account stellar winds and assume that stars above 25 M form 
3 M Q black holes after explosion, which have each a suitable lo- 
cal supply of gas such that accretion adds energy to the gas at a 
rate of 20 % of the Eddington luminosity. The stars between 10 
and 25 M Q are assumed to form 1 .5 M neutron stars, which also 
contribute 20 % of their Eddington luminosity. In a second sce- 
nario [no BH, late SN), we assume that stars with initial masses 
> 25 M B do not explo de, but directly form black holes (compare 
IDecressin et al.ll2()10i and references therein). The only energy 
sources are now SNe originating from stars with M < 25M Q . The 
third scenario (dark remnant accretion) assumes sudden accre- 
tion onto all black holes, which have formed as a result of the 
stellar evolution of the stars with M > 25M , accompanied by 
an energy transfer of 20 % of the Eddington luminosity to the 
gas, and no other energy source. Thus this scenario is applica- 
ble to the epoch when all core-collapse SNe have already taken 
place. As a subcase, we add accretion onto all neutron stars with 
the same efficiency. For these assumptions we integrate Eq. ([]]) 
using a fourth-order Runge-Kutta method with a sufficient time 
resolution to reach numerical convergence. 

2.2. Shell kinematics 

When neglecting gravity, the general analytic solution of the 
spherically symmetric thin-shell m odel is known for arbitrary 
mass profiles and energy input laws dKrausdl2003l) . With power 
laws for density (p oc r K ) and energy injection (E(t) oc t d ), the 
bubble expansion law is obtained as 



d+2 

r oc t 



(4) 



The gravitational pull peaks at r c / y2, and approaches zero for 
large radii. Since the density in the Plummer model drops like 
r , the exponent of t in Eq. © approaches infinity in the 
limit of large r for any reasonable (Fig. [TJ energy exponent 
d > -2. Hence, we generally expect that the superbubbles are 
fairly slow around n/2, where the gravitational pull is strongest, 
and quickly accelerate once they have overcome the gravita- 
tional potential well. This triggers the RT instab ility. Whenever 
a - g > 0, the instability grows on length scales (IChandrasekhaJ 
119611: iBernstein & Bookf l978l) A = (a-g)r 2 , where t is the time 
the instability is given to grow. Since g < 0, only a deceleration 
stronger than the gravitational acceleration may stabilise an ex- 
panding shell. This is t he case for many st andard situations of 
interstellar bubbles (e.g. lWeaver et alll977l) . Yet, acceleration is 
unavoidable in the case we consider here. In the following, we 
make the simple assumption that the bubble has burst and the 
pressurised hot gas escapes from the cluster, when A reaches the 
bubble radius. The shell material would then collapse back into 
the cluster, unless it has already reached escape speed. We take 
the time during which the bubble shows significant acceleration 
to define t. 

We show the energy injection together with the resulting 
bubble kinematics for the three respective assumptions about the 
energy injection law in Fig. Q] In the standard case (all energy 
sources active), the shell expands very slowly to about 4 Myr 
and in an oscillatory manner, as seen from the alternating sign of 
the acceleration. Once it reaches 7*1/2, the acceleration increases 
strongly and stays at a high level throughout, as expected from 
the above analysis. The shell reaches the local escape speed (red 
dashed line in the middle panels) at about 4.3 Myr, which is too 



2 



Martin Krause et al.: Superbubbles in globular cluster infancy I 



Cumulative energy output, initial mass: 9 Mio Msun 



Shell kinematics: 9 Mio Msun - standard 
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Shell kinematics: 9 Mio Msun - no BH acc, late SN 



Shell kinematics: 9 Mio Msun - dark remnant acc 
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Fig. 1. Produced energy for the standard scenario (top left) and superbubble kinematics for different assumptions about the energy contributors; top right (standard): 
winds and SNe for all massive stars and energy output from all black holes and neutron stars; bottom left (no BH acc, late SN): only SNe of stars less massive 
than 25 M e that explode after 8.79 Myr; bottom right (dark remnant acc): only sudden accretion onto the dark remnants (thick lines: black holes, only; thin lines: 
also neutron stars). The timescale for the global evolution of the GC is chosen at the birth of a coeval first generation of stars. The abscissae indicate the respective 
starting times of the three considered ejection scenarios, which corresponds to the moment when the considered energy sources become available. In scenario 3 this 
could happen any time once all the dark remnants have formed, i.e., after the last SN at 35 Myr after the birth of the 1G stars. Within each kinematics plot, the upper 
diagram shows the bubble radius (solid line) and the Rayleigh-Taylor scale (dash-dotted line), with the red dashed line indicating the half-mass radius. The middle 
diagram displays the shell velocity (solid line) and the escape velocity at the current bubble radius (red dashed line). The acceleration (positive: solid black line, 
negative: solid blue line) is shown in the lower diagram, with the gravitational acceleration at the current radius shown as a red dashed line. 



late to avoid pressure loss due to the RT instability. The cluster's 
gas, which is now in the shell fragments, should therefore re- 
main bound to the cluster. For the second scenario (no BH, late 
SN), the shell starts to fast accelerate after about 0.5 Myr after 
the SN activity is assumed to start. However, the escape speed is 
reached only after 1 .2 Myr. The shell is RT unstable long before, 
and therefore this scenario does not lead to gas expulsion, either. 
The situation is different for the third scenario, dark remnant ac- 
cretion: Here, the shell reaches escape speed immediately, due 
to the sudden power increase. The velocity drops slightly be- 
cause of the hydrodynamic evolution of the bubble. The escape 
speed is reached again after only 0.06 Myr and even 0.03 Myr 
if one includes the neutron stars. The RT instability is not able 
to affect the entire shell, and consequently, the gas is expelled 
from the cluster. We have also investigated cases with initial 
masses M tot = 10 6 M Q and M tot = 2 x 10 7 M Q with very similar 
results. The only difference is that for the high-mass cluster, gas 
expulsion by dark remnant accretion only works with the help 
of the neutron stars. The crossing time for the model clusters 
t c = 2.82r^2(GM tot ) 1/2 , is 0.22, 0.05, and 0.034 Myr, for the 1, 
9, and 20 xlO 6 M Q clusters, respectively, using the definition of 



iDecressin et al. For all our dark remnant accretion cases 

with the exception of the black-hole-only case for the high-mass 
cluster, the shells reach the half-mass radius much faster, and 
we may expect that the o uter 1G stars will also be lost (com- 
pare |Decressin!etjn]|20T(J. We have investigated the parameter 
space in e s f and r\ ji. For both parameters, we find critical values, 
above which gas expulsion by SNe only, or by the power sources 
in our standard scenario would be possible. They are excessive, 
apart from perhaps the lowest mass case, where a GC may lose 
its gas at e s f > 47 %. Yet, th is value is too high to expell a sig- 
nificant number of 1G stars ( IDecressin et al. l l20Tol) . Similarly at 
e s f = 0.33, the SNe succeed only for n/2 > 4 P c > which is on the 
high end of the observed values. 

3. Gas expulsion powered by dark-remnant 
accretion and possible observational tests 

Up to now, the gas expulsion scenario via SNe has been cen- 
tral to the two main scenarios for self-enrichm e nt in GCs (e.g. 
i Baumgardt et al.l 120081 : iD'Ercole et alJ 120081: IDecressin et al.l 
120101) . Here we show that this does not generally work for simple 
and standard assumptions about gas expulsion via a superbubble. 
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Table 1. Minimum star formation efficiency and half-mass ra- 
dius for SN feedback to be able to expel the gas 



M tot fl /10 6 W o 


e sf,cnt,std 


^sf,crit,sn 


n ji,crit d /pc 


1 


0.97 


0.47 


4 


9 


>0.99 


>0.99 


18 


20 


>0.99 


>0.99 


35 



Notes. Total initial gas mass out of which the cluster forms. 
^ Critical star formation efficiency for our standard energy injection 
scenario and a half-mass radius of 3 pc. ^ Critical star formation ef- 
ficiency for the 'SN only' energy injection scenario and a half-mass 
radius of 3 pc. Critical half-mass radius in the 'SN only' scenario, 
assuming a star formation efficiency of 0.33. 



While the energy injected by SNe in total is sufficient, it is not 
delivered fast enough to overcome the RT instability. The result 
should not be restricted to the Plummer model: Any GC forma- 
tion scenario should involve a strongly concentrated gas cloud. 
The gravitational pull will always be strongest on scales com- 
parable to the half-mass radius. Sudden acceleration, when the 
gravitational well is overcome, and RT instability, are the natural 
consequences. The asymptotic acceleration is particularly strong 
for the Plummer model. It is likely to occur, however, albeit at 
a weaker level for all reasonable profiles (compare Eq. [4] and 
Fig.ID. 

The only way to overcome the shell destruction by the RT 
instability is to inject the energy sufficiently fast, such that the 
gravitational well becomes less important for the dynamics. We 
show that a sudden activation of all dark remnants - a hypoth- 
esis, details of which need to be worked out in the future - is 
plausibly sufficient for this purpose. Here, we have adopted a 
general efficien cy factor of 20 %, ma inly for consistency with 
previous work dDecressin et al.ll2010l) . How may we motivate 
such an efficiency factor for dark remnant accretion? The Bondi 
accretion rate, lO~ 4 M^n^/v^M /yi, M3 being the remnant mass 

in units of 3M Q , «6 the ISM number density in 10 6 cm -3 and 
vi the higher of relative velocity and ISM sound speed in km/s, 
may exceed the Eddington accretion rate by a large factor when- 
ever the star is close to its outer turning point. Assembling mass 
to its vicinity in this phase, a remnant may plausibly be acti- 
vated for a large part of the orbit. Accretion onto compact stel- 
lar sized objects u sually leads to emi ssion in the X-ray part of 
the spectrum (e.g. IChiang et al.ll2.Q10b . For our standard model 
cluster, the hydrogen column density is about Nh ~ 10 25 cm -2 , 
i.e. it is Compton-thick and a large part of the radiated energy 
might be absorbed. The shells should ha ve densities of about 
10 6 " 7 cm -2 and accordingly cool down (ISutherland & Dopital 
119931) . collapse, remain neutral, and therefore also absorb X- 
rays efficiently. Additionally, there might be jets, which in the 
case of supermassive black holes are also know n to come some- 
times close to the Eddin gton luminosity (e.g. iRrausd l2005bt 
ICelotti & Ghisellinil2008l) . Since the compact objects in this sce- 
nario would not accrete from a binary companion, but from the 
general ISM, one may speculate that the jet powers might also be 
comparable to the supermassive black hole case. Jets communi- 
cate their energy efficien tly to their surroundings via radio lobes 
(e.g. iGaibler et al.ll2009h . " For our dark-remnant accretion sce- 
nario, we find a limiting initial cloud mass of about 10 7 M o above 
which the cold gas may not be ejected and therefore might form 
additional stars. This might contribute to explanations of the ob- 
served d ifferences at the high -mass end of GCs (e.g. greater Fe 
spreads, ICarretta et al.ll2010T) . 



It would be challenging to detect the dark remnants dur- 
ing their active phases, as the emission peaks in the X-ray part 
of the spectrum and the clusters should be Compton-thick at 
that time. The total X-ray luminosity of the cluster should be 
low, about 10 41 erg/s, the active time from our calculations 
is only about 10 4 - 10 5 years, and the prime objects of in- 
terest are high-redshift galaxies. If the radio luminosity were 
similar, one would expect fluxes of about fiJy, well in the 
reach of the upcoming Square K ilometre Array radio telescope 
(SKA) (e.g. iKrause et alJ I2009L and references therein). One 
could at least constrain well-defined models of cluster formation. 
If GC-formation was triggered by galactic scale shock waves 
(lHarris & Harris! 201 ll). e .g. associated with galactic winds and 
iets dKrausdl2002i ' 2005 ah . where the GCs are supposed to form 
in a galactic wind shell, one might expect some active GCs dur- 
ing the time when the jet is active too, thus marking the relevant 
evolutionary epoch of the GCs. The dark-remnant driven bub- 
bles could easily reach sizes of kpc or even 10 kpc on timescales 
of about 10 7 years. They might then leave a signature in the 
Faraday rotation signal if seen against a polarised background 
source. The SKA should also be able to dete ct high-redshift ra- 
dio lobes in polarisation (IKrause et al.ll2009h . The plasma clos- 
est to the radio sources is usually responsible for a big part 
of the rotation measure. Changing the structure of this ma- 
terial can leave observable features in Faraday rotation maps 
dHuarte-Espinosa et alj|201ll) . 
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